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Overall, brane inflation is compatible with the recent analysis of the WMAP data. Here we 
explore the constraints of WMAP and 2dFGRS data on the various brane inflationary scenarios. 
Brane inflation naturally ends with the production of cosmic strings, which may provide a way to 
distinguish these models observationally. We argue that currently available data cannot exclude a 
. . . non-negligible contribution from cosmic strings definitively. We perform a partial statistical analysis 

' of mixed models that include a sub-dominant contribution from cosmic strings. Although the data 

I favor models without cosmic strings, we conclude that they cannot definitively rule out a cosmic- 

. string-induced contribution of ~ f 0% to the observed temperature, polarization and galaxy density 

^ I fluctuations. These results imply that Gfi < 3.5 x 10"^(A/0.25)^/B/0.f , where A is a dimensionless 

5-H , parameter related to the interstring distance, and B measures the importance of perturbations 

induced by cosmic strings. We argue that, conservatively, the data available currently still permit 
■^L ' B OA. Precision measurements sensitive to the B-mode polarization produced by vector density 

perturbation modes driven by the string network could provide evidence for these models. Accurate 
determinations of ns{k), the scalar fluctuation index, could also distinguish among various brane 
inflation models. 
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PACS numbers: 98.80.Cq 



I. INTRODUCTION 



o , 

CO , Observations of the cosmic microwave background (CMB) P, J3 support the idea that the standard big bang 
' phase of the expansion of the Universe was preceded by inflation Recent results from the Wilkinson Microwave 
, Anisotropy Probe (WMAP) constrain the properties of proposed inflationary models tightly, but although 

■ some models are now excluded, numerous possibilities remain. A further challenge to observational cosmology is to 

CIhI try to hone in on a small class of viable models, even if identifying a single, correct theory of inflation may prove 

^ i impracticable. 

'"7^ ' All of the data collected up until now are consistent with a relatively pristine Universe in which the perturbations 
^ . observed today result from the amplification and distortion of a relatively featureless, Gaussian spectrum of fluctu- 
' ations produced by quantum effects during inflation. However, it is likely that inflation itself could have left behind 
, other remnants - such as cosmic strings - which could actively perturb both the CMB and dark matter of the Universe 
d • up to the present day. 

It is well-known that cosmic strings cannot be wholly responsible for either the CMB temperature fluctuations or 
the observed clustering of galaxies roughly speaking, the implied limits on the cosmic string tension fi allowed 
by observations is ^ 10~^. However, now that cosmology has entered an era in which the properties of the 
Universe are being revealed to unprecedented precision, a natural question is to what extent the observations can 
allow previously unwanted ingredients, such as cosmic strings (e.g. p^). Indeed, as the precision of cosmological 
observations increases, we might hope to be able to distinguish among numerous presently viable models for inflation 
by the properties of the cosmic strings they predict. 

Although the idea that inflationary cosmology might lead to cosmic string formation is not new 11], it has received 
new impetus from the brane world scenario suggested by superstring theory. In brane world cosmology, standard model 
particles and interactions correspond to open string (brane) modes, while the graviton, the dilaton and the radions 
are closed string (bulk) modes. Thus, our 3D Universe can be thought of as residing on a brane or stack of branes 
with three dimensions of cosmological size. These branes in turn reside in extra dimensions that are compactified. In 
such a model, inflation can result during the collisions of branes coalescing to form, ultimately, the brane on which 
we live [l^. 

In these brane inflation models, the separations between branes in the compactified dimensions are scalar fields 
(open string modes) that can act as inflatons, with the interaction potential between spatially separated branes 
providing the inflaton potential. Details of the brane inflation scenario depend on both qualitative and quantitative 
features, such as whether collisions involve a brane-antibrane pair f'^o branes coalescing at an angle 0|, as 
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well as parameters such as the sizes of the compactified dimensions |l5l. Il6|. Qualitatively, though, it appears easy 
to find models that predict adiabatic temperature and dark matter fluctuations capable of reproducing all currently 
available observations. A seemingly unavoidable outcome of brane inflation, though, is the production of a network 
of cosmic strings whose effects on cosmological observables ranges from negligible to substantial, depending 

on the specific brane inflationary scenario |18| . 

Although cosmic string production towards the end of inflation is possible in field theory models the scaling 
properties of the cosmic string networks in brane inflationary scenarios are different than that in the familiar 3+1 D 
simulations, since intercommutation probabilities are smaller as a consequence of the existence of extra dimensions. 
In addition to placing constraints on the amplitude of string induced perturbations of the CMB, we show that the 
results place limits on G^/X, where n is the string tension and A < 1 is related to the interstring distance as a fraction 
of the horizon. The value of A is sensitive to the intercommutation rate of strings, however the exact relation is 
expected to be model-dependent. 

Here we shall first review the essential points of brane inflation, and examine the constraints imposed by the 
WMAP observations if we ignore the contribution from cosmic strings. These constraints allow us to delineate a 
range of possible cosmic string tensions. Then, we assess quantitatively the extent to which cosmic strings can 
contribute to the CMB temperature fluctuations and power spectra of dark matter density perturbations. In this 
analysis, we hold properties of the background cosmological model fixed to their best fit values, as determined by 
WMAP "5^ without cosmic strings. (A more detailed analysis that varies the background cosmology as well is under 
way.) Although the available data favor models without cosmic strings, they may still allow, within the uncertainties, 
a contribution from string- induced perturbations of up to 10%. They also imply scalar perturbation indices ns{k) 
which, although still consistent with a broad range of models, may be able to discriminate among them in future. We 
also compute the dark matter density perturbation power spectrum, and compare with observational determinations 
from the 2dFGRS galaxy survey [13 . We discuss the interpretation of these results in terms of the string tension 
and efficiency with which the cosmic string network decays via intercommutation of string segments, which is reduced 
in a Universe with extra dimensions. Finally, we discuss the prospects for detecting B-mode polarization, which is 
expected to be a prominent signature of a cosmic string network, in view of the constraints implied by our analysis. 



II. BRANE INFLATION AND COSMIC STRING PROPERTIES 



Recently, the brane world scenario suggested by superstring theory was proposed, where the standard model of 
the strong and electroweak interactions are open string (brane) modes while the graviton and the radions are closed 
string (bulk) modes. The relative brane positions (i.e., brane separation) in the compactified dimensions are scalar 
fields that have just the right properties to act as inflatons. Thus, the brane inflation scenario emerges naturally in 
the brane world In this picture, the inflaton potential is due to the exchange of closed string modes between 

branes; this is the dual of the one-loop partition function of the open string spectrum, a property well-studied in string 
theory po| . This interaction is of gravitational strength, resulting in a very weak (flat) potential, ideally tailored for 
inflation. 

The potential is essentially dictated by the attractive gravitational (and the Ramond-Ramond) interaction between 
branes. As the branes move towards each other, slow-roll exponential inflation takes place. This yields an almost 
scale- invariant power spectrum for density perturbation, except there is a slight red-tilt (at a few percent level). As 
they reach a distance around the string scale, the inflaton potential becomes quite steep so that the slow-roll condition 
breaks down. At around the same time, a complex tachyon appears, so inflation ends rapidly as the tachyon rolls 
down its potential. In effect, inflation ends when the branes collide, heating the universe to start the standard big 
bang phase of cosmological expansion. This brane inflationary scenario may be realized in a number of ways p^ll6|. 
The scenario is simplest when the radion and the dilaton (bulk) modes are assumed to be stabilized by some unknown 
non-perturbative bulk dynamics at the onset of inflation. Since the inflaton is a brane mode, and the inflaton potential 
is dictated by the brane mode spectrum, it is reasonable to assume that the inflaton potential is insensitive to the 
details of the bulk dynamics. 

Coupling of the tachyon to inflaton and standard model fields can allow efficient heating of the universe if certain 
conditions on the coupling of the tachyon to standard model particles are met ^^l^ . As the tachyon rolls down its 
potential, besides heating the universe, the vacuum energy also goes to the production of defects, in particular, cosmic 
strings. The effect of the resulting cosmic string network may be negligible or rather substantial, depending on the 
particular brane infiationary scenario |l8j . However, in all cases, we expect the density perturbation power spectrum 
in CMB to be dominated by the adiabatic fluctuations arising from quantum fluctuations of the inflaton during 
brane inflation, not by the nonadiabatic contributions from cosmic strings. However, the contribution to the density 
perturbation power spectrum in CMB coming from the cosmic string network may be large enough to be observable. 

We devote this section to a review of the implications of brane inflation. For a broad set of models, we present 
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results for the slow roll evolution, fluctuation spectra, string mass scale, and associated cosmic string tension. (These 
results follow directly from the treatments in Refs. [iMliailiflllSilill-) We consider the coUision of a Dp brane with 
a Dp brane at an angle 9\ collision with a Dp antibrane corresponds to 9 — t:. (Here and throughout this section, we 
follow pi^, which contains more details and discussion). Of the ten spacetime dimensions, one is the time, three are 
the large spatial dimensions we live in, and the rest are compactificd. Of the compact dimensions, p — 3 are parallel 
to the brane, and we take their compactification lengths to be l\\ = 27rr||, implying a volume Vj| — ^jj"^- Of the 
remaining d = 9 — p dimensions, we take d — d± to be compactified with a size 2?:/ Ms, where Ms is the string scale, 
while the remaining dj_ are compactified with a size £j_ = 2Trr± > 2'k/Ms. The 10-dimensional gravitational coupling 
constant is 

where gs is the expectation value of the dilatonic string coupling, which is related to the standard model gauge 
coupling a{r\\) on a scale l/r|| by 

= 2(M,r||)P-3a(r||) ; (2) 

the 4-dimensional Planck scale Mp — (SttG)^^/^ is then 

2.,2 _ M^s{Msr^Y-{Msr\\Y-^ 
gsiyip — ■ \o) 

TT 

The outcome of brane inflation will therefore depend on several parameters, p, d^, r^, r\\ and a(r||). 

We will distinguish between two different potentials for the interaction between branes, depending on their sep- 
arations. (See Hi El Hi.) For some scenarios, a fixed lattice of branes is considered to be spread throughout 
the compactified dimensions, with a moving brane placed inside one lattice square. At separations small compared 
to the lattice size of the compactification topology, the interaction is "Coulombic", with a potential of the form 
V{y) = Vo — U jy"^^^"^ for a separation y in the large compact dimensions. This potential is suitable for inflation 
resulting from the collision of a pair of relatively nearby branes at a small angle [lj| • When the separation is nearly 
equal to the lattice size, an expansion about zero displacement from the anti-podal point gives T^(j/) = Vb — U'y'' , 
where a depends on the compactification topology. This potential is suitable for the brane-anti-brane scenario (which 
corresponds to branes at angle tt). In the next two subsections, we summarize the infiation scenario for interbrane 
potentials of these two general forms. 



A. Coulombic Inflation 



Consider a potential of the form 



{di_ - 2)V"i^ 



(4) 



with i/) cx y, the interbrane spacing; for the special case d± = 2 this becomes a logarithmic potential, but the results 
we derive below may be applied to this special case. (We only consider dj_ — 2 > here to simplify our analysis, since 
the results generalize easily to the logarithmic case.) In the slow roll approximation, the equation of motion for ^ 
becomes 



dL 



■qMl 



(5) 



where i = Ina is the logarithm of the scale factor a(t), which we take to be zero at the start of inflation. The slow 
roll solution is then 



- d^T^MlL = [rf^r/A/l (Li„f _ L)] = {d^i^MlL,) 



where the starting value of the field is i/'i, the total number of e- folds in inflation, is 



d^rjMl ' 



(6) 



(7) 
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and Lr = L\, 



L is the total number of e-folds remaining in inflation. The curvature fluctuation spectrum is then 



ij4 



V^id^Lr) 



(8) 



where Lj. is evahiated when k/a = H, or ln(fc/fco) — ir,o ~ L^, where fcg is a reference scale, which crosses with L^^o 
e-folds remaining in inflation. The fluctuation spectrum is very flat, with only slowly varying scalar index, ng{k): 



ns{k) - 1 

dns{k) 
din fc 



dlnA2j(A:) 
din fc 
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Lr{k) 



-6 X 10" 



-0.03 



60 



60 



Lr{k) 



Lr{k) 



(9) 



both of which are in the range of uncertainty of the determinations in Q . 

The challenge to this, or any other, inflation model is to have sufficient inflation as well as small curvature fluctuation. 
Since the precise value of Lmf depends on initial conditions as well as parameters of the model, let us first consider 
the constraints on the latter implied by comparing Eq. (jHJl to the WMAP result A^(A:o) = 2.95 x 10"^A(/co) with 
^(^o) = 0.9 ± 0.1. To do this, let us consider a particular model with p = A and a small collision angle, 9\ then we 
have 



Td\\ 

2 



^ 2 

4 

327r'^Q;(r| 
167r 



M. 



6-d I 



{ed±Lr) 



24[64/3(di)]i tt"^^^! ["(^-Ii)]""'^ 



Mp 



where 



1 1 ^ . _ o 



Let us consider the specific example d^ — 2; for this case we find 

4 



9Lr 



A/, 



7687r6[a(r||)]2 \Mp 



[d± = 2] , 



and therefore the string scale is determined to be 



Wp 



^ ^ 2.5 X 10-2[25a(r||)]V2[A(fc„)]i/4 ( _^ 



10 



1/4 



[di_ = 2] 



that is, on the same order of energy as the GUT scale, 10 GeV. Larger d± leads to smaller Mg/Mp; thus if d± 
we find 



Alik) 



{0Lr) 



3/2 



12^/27^7/2a(r||) \Mp 



Ms 



[d± = 4] 



(10) 

(11) 
(12) 

(13) 

(14) 
4 

(15) 
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which in turn requires 



^ c.1.6 X 10-3[25a(r||)]i/3[A(fco)]i/=^ ( ^ 



1/2 



The total number of e-folds in inflation is 



±(d±-l) 

d I +2 ( 1/-. A,f ^ ^a 



iinf — 



647r5/3(di)6l[a(r||)]2Af|, 



3eij^ 2(dj_-l) 



0.025[A(fco)]i/2 / 10 



1/2 



[25a(r||)](106') V^*^-- 
0.025[v4(fco)]^/3 /^O 
[25a(r||)]4/3(l06l) \eL, 



^ ] [QMsT^f = 2] 



[d^ = 4] 



3A^(fc) 



8/3(dj.)a2(n| 



(16) 



(17) 



where we have let yi — 27Tr±Q with d < 1. To get Li^i ^ 60, we must require C,iMsr± > 50 for dj^ = 2 or CiMsr± ^ 10 
for d± = 4. Note, though, that for large 9, it is not possible to have enough expansion during inflation. In this case, 
the images of one brane exert non-trivial forces on the other brane, resulting in a power-law type potential. 



B. Power Law Inflation 



Next, we consider potentials of the form 



(18) 



such potentials arise for a brane situated near the origin. The value of a depends on the compactification topology. 
For hypercubic compactification, = 4, whereas in other cases, a = 2. Note that in actuality the potential need not 
depend just on interbrane separation in such a picture, and the trajectory of the brane can be complicated. Here, 
though, we confine ourselves to simple one dimensional (diagonal) brane motion. 

Following Eq. IjlSI) . we see that the origin - = - is an unstable equilibrium point, and any perturbation away 
from it will result in slow motion of the brane. For a > 2, the slow roll solution is 



1p = [■0, 



(T-2 



a{a - 2)r/Mj,L]— [cr(cr - 2)r]Mj,Lr]~ , 



and the total number of e-folds in inflation is 



Jinf 



(19) 



(20) 



where is the starting value for the inflaton. Quantum fluctuations will imply = Ci_ff/27r, where Ci ^ 1. The 
curvature fluctuation spectrum is 



A^(fc) = 



The implied fluctuation spectrum is acceptably flat: 



Vo{(J - 2f[a{a - 2)t]] — M 



2(4-g) 2(<T-1) 

cr-2 r CT-2) 
P J^7- 



127r2 



71s(fc) - 1 = 



dlnA^(fc) 
dink 

dns{k) 
din fc 

For fJ = 4, Eqs. if^ and ((201 become 



2(^-1) 
{a~2)Lr{k) 

2(^-1) 
Ll{k){a-2) 



0.03(ct- 1) 



cr-2 

-6 X 10""* 



60 



Lr{k) 
(7-1 



cr-2 



60 



Lr{k) 



A^(fc) = 
Lint = 



CfVo 
96tt2MI,{t]Vo) 



(21) 



(22) 



(23) 
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the observational constraints on the curvature fluctuation spectrum therefore require 

Stt^A^ (fc) , /60\^/' 



vVo ^ _;}^ ^ - 2.5 X lO-^M(fco) — 1 , (24) 



i.e. the potential must be extremely flat. This requirement is well-known from studies of new inflation, which 
sometimes idealize the potential to Eq. p8|l with a small dimensionless parameter A equivalent to 77 Vq. In Ref. (161] . 
a particular toroidal compactification is proposed where this small parameter is {F is a geometrical factor related to 
the compactification geometry) 

^^0 - , (25) 



which can be small enough for 9 ^ 0.1 provided that 

In this picture, the flatness of the effective potential is attributed to a relatively small value of the string scale compared 
with the Planck mass. 

Special treatment is required for a — 2, which is expected for any non-hypercubic compactification topology. For 
this case, the scale factor grows like a powerlaw in time during slow roll: 

where V'/ and a/ are the values of the field and scale factor at the end of slow rolling. Since d\ntp/d\na = 2riMp, 
we require rjMp <C 1 for slow rolling. It is easy to see that for this potential, i/j/SHip = 2rjMp/3 ^ 1. Slow roll 
ends, for this potential, only when ?/'^/2Vb — 2rfAIpip'^/3 ~ 1, or ipf ~ {riMp)~^ ^ Mp, or when the polynomial 
approximation to the potential fails, which happens when the brane moves a substantial fraction of a lattice spacings. 
The total number of e-folds in inflation is 



^'"f - 2^Ml - 2^Ml ■ ^^^^ 



The curvature fluctuation spectrum for this case is 



2 _ {H/2^ijffe^^^'''^^^ ^ (g/2^V>,)^(a,/a)4'>^4 ^ {^H /2^ij.,f{a,H /kf^^l 

' ^{iiMlY ^{riMlY 4(r;M2)2 ' ^^•^> 



where evaluating at horizon crossing implies that k — Ha, which has been used to get the final form of the spectrum. 
In this case, 

dlnA^(fc) o , , 

ns-l= ^^^^ ' = -4?7M| , (30) 

which is independent of k. The WMAP analysis implies that yyAfJ, ^ 0.01. From the first form of Eq. (|29|l . and 
■qAIp ~ 0.01, it follows that the observed temperature fiuctuations can be accounted for if (Vo = 2t4£||) 



27rV'/ ' V 0.01 

in which case Linf ^ 10'^. For yf ~ 2TTr±, this relation implies 



(2^Mp2;^)-^~10-M , (31) 



Mp V 0-01 

which is generally smaller than our previous estimates unless M^r^ ^ 10^ 



- 10-^ ( ^ ) Msr. , (32) 
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C. Cosmic String Properties 

Because the inflaton and the ground state open string modes responsible for defect formation are different, and 
the ground state open string modes become tachyonic and develop vacuum expectation values only towards the 
end of the inflationary epoch, various types of defects (lower-dimensional branes) may be formed. A priori, defect 
production after inflation may be a serious problem. Fortunately, it is argued in Ref . \wL Il7j that, from the properties 
of superstring/brane theory and the cosmological evolution of the universe, the only defects copiously produced are 
cosmic strings. In superstring theory, Dp-branes come with either odd p (in IIB theory) or even p (in IIA theory). 
The collision of a Dp-brane with another Dp-brane at an angle (or with an anti- Dp-brane) yields D{p — 2)-solitons 
(i.e., codimension 2). Topologically, a variety of defects may be produced. Because they have even codimensions with 
respect to the branes that collide, they have specific properties |22l |. Cosmologically, since the compactified dimensions 
tangent to the brane are smaller than the Hubble size, the Kibble mechanism works only if all the codimensions are 
tangent to the uncompactified dimensions. As a consequence, only cosmic strings may be copiously produced [T^ll7j|. 

The observational imprint of cosmic strings is determined primarily by the product of Newton's constant and the 
cosmic string tension, Gfi, assuming the evolution of the string network can reach the scaling regime. The value of n 
implied by superstring cosmology depends on several parameters, but is most sensitive to the string scale, M^. To get 
an order of magnitude estimate, we may use the small case, which is arguably the most likely inflationary scenario. 

The cosmic strings may be Dl-branes, but most likely, they are D{p — 2)-branes wrapping around (p — 3)-cycles in 
the compactified dimensions. If the Dl-brane is the cosmic string (i.e., p = 3), its tension is simply the cosmic string 
tension: 

M = n = -^. (33) 
2ngs 

However, we expect the string coupling generically to be > 1. (It is well-known that radion and dilaton moduli are 
not stabilized by perturbative dynamics in string theory. Presumably, any superstrongly coupled string model is dual 
to a weakly coupled one, and thus cannot stabilize the moduli either. We therefore expect a moderately strong string 
coupling, since only in this case will we find non-trivial dynamics.) To obtain a theory with a weakly coupled sector in 
the low energy effective field theory (i.e., the standard model of strong and electroweak interactions with weak gauge 
coupling constant a), it then seems necessary to have the brane world picture, in which we have the Dp-branes for 
p > 3, where the (p— 3) dimensions are compactified to volume V\\. Now the cosmic strings are D(p— 2)-branes, with 
the {p— 3) dimensions compactified to the same volume Vj[. Noting that a Dp-brane has tension Tp = M^^^ /{27T)Pgs, 
the tension of such cosmic strings is 

^ {2TT)P-^g, 2'Kgs 4a7r " ^ ' 

for a ~ OiGUT — 1/25. For one pair of branes at angle 0, only this type of cosmic strings is produced topologically. For 
a large enough stack of branes colliding, the Dl-branes may also be allowed topologically, but they are not produced 
cosmologically. Thus, /i ~ 2M^ is a reasonably general estimate. We considered estimates of Ms implied in various 
scenarios for brane inflation in S III Al and IIIBI These estimates are broadly consistent with 

IQ-^ >Gfi> 10"" , (35) 

although a smaller range is obtained in any specific model or class of setups. For example, for branes colliding at a 
small angle, a likely range is 

5 X lO"'^ > > 7 X 10"^ . (36) 
Thus, brane inflation can lead to cosmic string tensions below, but not far below, current observational bounds. 

D. Tensor modes 



During slow roll, the tensor power is 



128G^K. 2Vo 

^'^^''>- 3 ~37r2M4 ^^^^ 
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and is smaller than the scalar power by the factor 



V 



M|, Vdlna 



(38) 



How small r{k) is depends on the specific brane inflation model. For branes intersecting at an angle 9 we find that 

2Vb O^M^ 



3n^Ml, 967r5a(r||)Mj, 
3.3 X IQ-^^ (109)^ [25a(r||)] 



lOA(fco) 

6Lr 



[d± = 2] 



5.4 X lO-i^(lO0)2[25a(r||)]i/3[v4(fco)]^/^ 
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6Lr 



[d± = 4] 



(39) 



In this case, the amplitude of the scalar mode is smaller than the amplitude of the perturbations due to cosmic strings 
by a small numerical factor times 9^, unless cosmic string intercommutation is extremely inefficient; see S III Cl and lllll 
below. For powerlaw brane-antibrane inflation, 6* tt, and Vo — M^/(27r)'^a [T3llTfil |. so for this case wc find 



127r5aM|, 



(40) 



Nominally, these perturbations can be comparable to those induced by cosmic strings, although they may be relatively 
suppressed by the small numerical factor (127r^Q;)"^ ~ 0.007/(25a). However, the spectrum of fluctuations produced 
by cosmic strings will still distinguish them from those due to primordial tensor modes. Both strings and the primordial 
tensor modes result in the B-type polarization of the CMBR. The predicted angular power spectrum, C^^ , has been 
calculated for tensor modes from inflation (see e .g. H^). It has a generic feature that most of the power is on larger 
angular scales, in the region / ^ 100. This is very different from the shape of the Cj^^ spectrum predicted by cosmic 
strings. There the dominant contribution comes from the vector modes and, as one can see from Fig. 01 most of the 
power is on smaller scales: 700 < I ^ 1000. 

As of today, the B-type polarization has not been detected |3 and the experimental constraint on r{k) is rather 
mild: r{ko = 0.002Mpc-^) < 0.71 ^. 



III. CMB, MATTER DENSITY AND POLARIZATION POWER SPECTRA: CALCULATION AND 

INTERPRETATION 

The fluctuations expected to arise from brane inflation should be an incoherent superposition of contributions from 
adiabatic perturbations initiated by curvature fluctuations A^(fc) and active perturbations induced by the decaying 
cosmic string network. For example, the resulting CMB temperature maps will yield 

Ci = VFCf'^'^''^"'^ + BCf ""'5", (41) 

were W and B are weighting factors. Analogous expressions hold for matter density and polarization power spectra. 
In Eq. (gU, the weight factors W and B determine the relative importance of the adiabatic and cosmic string 
contributions. We choose the weight factor W so that W = 1 when there are no cosmic strings. 

In computing the combined effects of adiabatic and cosmic string perturbations, we have kept the cosmological 
background parameters flxed at their best-flt values according to j^J. In addition to B and we vary the spectral 
index of the scalar curvature fluctuations, Us. The tensor contribution to the adiabatic component was set to zero, 
since, as discussed in § III Dl it is likely to be small. When fltting to both WMAP and the 2dFGRS data, we considered 
two cases (described in more detail in S lIIIBjl : with bias b fixed, and with b being an additional parameter of the fit. 



A. Cosmic strings and the CMBR 

Perturbations due to cosmic strings were calculated using the model first introduced in and further developed 
in ll^ll^. The main idea is to represent the cosmic string network by a collection of uncorrelated, straight string 
segments moving with random, uncorrelated velocities. All segments are produced at some early epoch and, at every 
subsequent epoch, a certain fraction of the number of segments decays in a way that maintains network scaling. The 
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length of each segment at any time is taken to be equal to the correlation length of the network which, together with 
the root mean square velocity of segments, are computed from the velocity-dependent one-scale model of Martins and 
Shellard |27j |. The positions of segments are drawn from a uniform distribution in space and their orientations are 
chosen from a uniform distribution on a two-sphere. 

This model is a rather crude approximation of a realistic string network. However, with a suitable choice of model 
parameters, its main predictions for CMB and matter power spectra have been shown to be in agreement with 
results obtained using other local string sources p^ . 129^ . The main advantage of our model is its flexibility. For 
example, parameters can be chosen to describe stings with different scaling properties, different amounts of small 
scale structure etc. This is especially valuable when describing strings produced in brane inflation, since strings are 
expected to intercommute with a lower probability in the presence of extra spatial dimensions |18| |. 

It is well known that properties and possible observational signatures of global and local strings can be dramatically 
different. Global strings predict almost no power on small angular scales for CMB temperature anisotrop y ISOl . while 
local strings produce a quite significant broad peak at I ~ 450 in a spatially flat universe |^ |28L 1291 l3lj . Also, 
local strings develop a signification amount of small-scale structure over the course of their evolution, which tends to 
suppress the vector component of metric perturbations. For that reason one would expect the strength of the _B-type 
polarization [s^ sourced by g lobal strings to be somewhat larger than that of local strings [s^l ■ 

Numerical simulations |34l l35l Is^ show that during the radiation and matter dominated eras the string network 
evolves according to a scaling solution, which on sufficiently large scales can be described by two length scales. The 
first scale, £,{t), is the coherence length of strings, i. e. the distance beyond which directions along the string are 
uncorrelated. The second scale, L{t), is the average inter-string separation. Scaling implies that both length scales 
grow in proportion to the horizon. Simulations indicate that ^(i) ^ t, while L{t) = jt, with 7 « 0.8 in the matter 
era 0, [S^ • The one-scale model |3^ , in which the two length scales are assumed to be equal, has been quite 
successful in describing the general properties of cosmic string networks inferred from numerical simulations. These 
simulations have assumed that cosmic strings would reconnect on every intersection. It is of interest to us, however, 
to also consider the case when the reconnection probability is less than one. If strings can move and interact in extra 
dimensions then, while appearing to intersect in our three dimensions, they may actually miss each other. Hence, the 
effective intercommutation rate of these strings will generally be smaller than unity. As a consequence, one would 
expect more strings per horizon in these theories. Because of the straightening of wiggles on sub- horizon scales due to 
the expansion of the universe, one would still expect S.{t) ~ t. However, the string density would increase, therefore 
reducing the inter-string distance. Hence, smaller inter-commutation probabilities imply smaller 7. 

Consider a string network in a volume V described by the two length scales introduced above. On average there is 
one string segment of length ^ per volume ^L^. The rough number of such string segments is 

N^^. (42) 

If the energy per unit length is /x, then the total energy of the string network is 

E = N^^^=^. (43) 

and the energy density is just p = E/V = f-i/L"^- 

Now suppose we want to calculate the effect of this string network on the CMB temperature anisotropy. In 
particular, we want to find the power spectrum, i. e. the 2-point function. For simplicity, let us assume that only 
p{t) affects CMB (in general we would have to consider all components of the string network's energy-momentum 
tensor T^). Then to evaluate the CMB power spectrum it suffices to know the 2-point unequal time correlators 
{p{k,ti)p{k,t2)) at all k, ti and t2- Here /9(k, t) is the Fourier transform of p(x, t). The CMB power spectrum is 
roughly given by 

Ci = j dk j dti J dt2 Li{h,t2,k)[{p{k,ti)pik,t2))] , (44) 

where I//(ti, t2, k) is a linear operator. 

Again, for simplicity, let us assume that f (i) ^ t and i = 7< at all times, namely that the network scales perfectly 
with time. We want to see how C/, or equivalently {p{k, ti)p{k, 12)), depends on 7 and /i. Let us also assume that the 
segments are straight. 

At time t there are roughly N[t) — Vji^L? segments and 

N{t) N{t) 

pik, t) = J2 P^^^ (k, t) ^ E ^^'^ *)^^(^) ' (45) 

1=1 i=l 
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where p^''\k,t) = p^'^\k,t) / {^£^{t)) to factor out dependences on /i and (^(t) (the phase of p may still depend on ^(t) 
but the amplitude does not). Now we can write: 

N(ti)N(t2) 

{p{kM)p{kM) ^ E <^*''(^'^i)'^^''(^'^2))M'e(ii)C(i2). (46) 

Because individual segments are statistically independent 

(pW(fc,ii)p(^)(fc,t2)) = 5.,j (p«(fc,ii)p«(fc,t2)) , (47) 

and therefore 

min[Ar(ti),Af(t2)] 

{p{kM)p{kM)) = E (P^')(fc,ii)/5«(fc,t2))/i2e(ii)e(i2) (48) 

i=l 

To interpret min[iV(ti), N{t2)] it might help to think that all segments were there at the initial time but over course 
of their evolution some of them decayed. For certainty, let us assume ti < t2 and therefore N{ti) > N{t2)- We can 
now write 

N{t2) 

{p{k,t,)p{k,t2)) = E (p^^Hfc,ii)p*^'(fc,<2))A^'C(ii)e(i2) =iV(t2)(p('nfc,^i)/5''Hfc,i2))/i'e(ii)^(i2) , (49) 

i=l 

where the last step is possible because all segments are statistically identical and hence 

(p«(fc,ii)pW(fc,t2)) = (/5(i)(fc,ti)p«(fc,i2)) foralH. (50) 
Substituting i^t,L = -ft and N{t) = V/{e,L'^) - V/{-iH^) we find 



V u? 

{p{k,h)p{k,t2)) = -^p^tit2{p^^\k,ti)p'^^\k,t2)) = ^ X F{k,tiM) , (51) 

7 ^2 7 



where F{k,ti,t2) is independent of ^ or p, and therefore 

^strings ^ P^ ^52) 

ry2 

Rather than working with the parameter 7 — L{t)/t, we define a related parameter A = 1{t)/t, where 1{t) is the 
comoving interstring distance, L/a, and r is the conformal time, dr — dt/a. In our code we take parameter A to be a 
constant. Correspondingly, the meaning of the parameter B is 



G(A,G/i) 



CiiXo,Gpo) 



pXo 
Mo A 



2 



(53) 



with Gpo = 1.1 X 10~^ and Ao = 0.25 adopted as reference values. 

The relation between parameter A and the string intercommutation probability p has been a subject of extensive 
research. Qualitative arguments in ;39] suggest that A cx p^/^. However, these arguments do not take into account 
possible effects of small-scale structure on strings. Results of the most recent numerical simulations suggest that 
for p > 0.1 the scaling string density is effectively the same as for p — 1, while for smaller p one should expect 
A oc p^'^. Such weak dependence of A on p can be explained by multiple opportunities for string reconnection during 
one crossing time, due to small-scale wiggles. 



B. Methods 



We have performed a partial statistical analysis in which we held the parameters of the background cosmological 
model (total, matter, baryon and dark energy density parameters, Hubble constant, reionization optical depth) 
fixed at their WMAP best fit values according to 0. More specifically, we considered a flat ACDM universe with 
^^CDM — 0.225, ilb — 0.045, JIa = 0.73, Hq = 71 km/s/Mpc and r = 0.17. The scalar spectral index Ug was allowed 
to vary within bounds set by the prior 0.8 < < 1.2, in increments of Aug = 1.25 x 10"'^. 
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The string spectra were calculated only once, using the string model parameters chosen to produce spectra that 
roughly agree with [S^. In particular, we set GfiQ = 1.1 x 10~^, which, if strings were the only source 

of inhomogeneity, would result in temperature anisotropy in a rough agreement with observations on COBE scales. 
The string intercommutation probability A was approximately A ^ 0.25 and was allowed to vary only insignificantly 
during the radiation-matter domination transition. 

The CMBR and linear matter power spectra for both adiabatic and string parts were computed using respectively 
modified versions of CMBFAST However, one cannot directly compare linear matter spectra outputted by 

CMBFAST to the galaxy clustering data published by the 2dFGRS team. We have "processed" the theoretical power 
spectrum P*'*(fc) for both adiabatic and string components following a procedure similar to that prescribed in Section 
5.1.4 of First of all, we output P*^{k) at the effective redshift of the 2dFGRS survey: z^s = 0.17 (the valued 
suggested in Q). Then we correct for the redshift space distortions using the approximate formula given in 

P^-^k) ^ Pf (fc) = {1 + ^3+ \l3'')P"\k) , (54) 

with (3 = 0.45. We then convolve it with the 2dF window function using the matrix Mw provided on the 2dFGRS 
website: 

Pf (fc) ^ Pw{k) = 9)^^(9) (55) 

and, finally, multiply by the bias factor to obtain the spectrum that can be compared to data: 

P(fc) = b^Pw{k) (56) 

where b = 

We have chosen to fit to the binned version of WMAP data, for both temperature (TT) and cross temperature- 
polarization (TE) angular spectra. When fitting to WMAP4-2dF we had a choice of making the bias b an additional 
parameter or using a prescribed fixed value. We have investigated both possibilities and will refer to them the "6-free" 
and "6-fixed" model. 



C. Results 



We are interested in constraining the fractional contribution to the WMAP and 2dFGRS data given by Eq. (|41|) . To 
do this properly, from first principles, would require using theoretical models with both cosmic strings and adiabatic 
perturbations to generate the relevant sky maps, and then compare these directly to the data to compute likelihood 
functions that can be used, along with appropriate priors, the posterior probability distributions for the parameters of 
the models. To get a quick and dirty bound on the most interesting parameter, B, we have chosen to adopt a simpler 
approach in which we treat the published results for C; of the TT and TE sky maps, and the 2dFGRS power spectra 
as data, and use them and their uncertainties to construct a statistic in the usual way. For each value of B, we 
minimize the statistic with respect to variations in the parameters n^, W, and the bias parameter, to find Xmini^)- 

If our "data" really represented independent data points, then we could be confident that at each B the value of 
Xmini^) would obcy a distribution that roughly has a mean value equal to the number of degrees of freedom in the 
fit, I/, and a standard deviation -\/2i/. Fig. Qshows the results for 5-free (left) and 6-fixed (right) (as defined at the end 
of S lmBll . For our procedure, we have N = 38(TT + 26(TE) + 32(2dF) = 96 "data" points. We keep cosmological 
paramters such as $7™, f^b, ^^A, h and r(reionization) fixed at their best fit values from Therefore, we have either 
three (n^, W ^ and b for 6- free evaluation) or two {ug and W for 5-fixed evaluation) parameters in each of the panels 
in Fig. n so J/ = 93 (5-freee) or 94 (6-fixed), and ^/iv ~ 14 in either case. Since the two panels in the figure actually 
show the reduced statistic 'x^{B)/v, we could regard values of B with 'x^{B)/v within approximately 0.15n of the 
minimum value as representing more or less equally good fits to the data to within n "sigma" . Below, we shall use 

/v ^ 1.25 for placing a limit on P; this can be interpreted as either a 1.7"cr" bound, or as a crude attempt to 
account for the fact that we are also ignoring the effects of varying several more cosmological parameters. Although 
assessing the acceptability of models with different values of B is not justified rigorously, we note that the minimum 
values of shown in cither panel of Fig.^are close to one, which may be taken as a posteriori assurance that our 
method for comparing models with different values of B may not be far off. 

The left panel of Fig. ^ corresponding to the 6-free model, contains plots of the minimized reduced Xxt/^i 
Xtt+te/^ and XTT+TE+2dFa/'^ ~ ^'^) ^ functions of B computed using, respectively, WMAP's TT spectrum 
only, TT and TE, and TT, TE and the 2dF spectrum. As one can see from that plot, in all cases the lowest value of 
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FIG. 1: (Left) Minimum reduced x^/i" 

as a function of B for the fit of the 6-free model to the WMAP's TT (red triangles), 
TT+TE (red stars) and TT+TE+2dF (red hollow squares) with corresponding best fit values of — 1 (green squares) and 
Us — 1 (blue circles). We only show W and Us from the fit to the TT data alone. The error bars correspond to the diagonal 
elements of the covariance matrix. Adding the remaining data sets changes the best fit values of W and ris by only a small 
amount, well within the plotted error bar. (Right) Same as on the left but for the fe-fixed model with all data sets included. In 
the left panel, u = 93 and in the right i/ — 94. The values of the reduced x^/*^ at B = are, in the left panel, 0.97 (TT), 1.1 
(TT+TE) and 0.98 (TT+TE+2dF), and 1.1 in the right panel. 



the minimum 'x^{B)lv occurs at i? = 0, i. e. when there is no contribution from strings. However, as additional data 
sets are added to the fit, the values of the minimum x^{B) become smaller. Fig. Q] also shows the values of Us — 1 
and W —1 values that minimized x^/'^ each B, with estimated uncertainties, of order ±0.01 for both. Not plotted 
are the values of the bias parameter b, which varied from 0.97 at B = to 1.02 aX B — 0.25. Note that the value of 
TT-s — 1 is typically around —0.02 ± 0.01 for B < 0.1, the region where the minimum x^{B)/v is less than about 1.25. 

The right panel of Fig. ^ corresponds to the 6- fixed model, i. e. the model in which the value of the bias was held 
fixed at & = n^/p. Results are shown for the minimum value of ~x^{B)/v constructed using all three data sets. Here 
we find that the lowest value of the minimum y^(B)lv is at i? = 0.025, not i? = 0. In view of the limitations of 
our statistical analysis, one cannot attach much significance to either the existence of this minimum, or the precise 
nonzero B at which it occurs, but it is worth observing that adding a small string contribution may be welcome when 
matching the CMBR normalized matter spectra to the galaxy data. Moreover, this result shows that at the very least 
we cannot exclude values of -B ~ 0.025. We note that the values of — 1 corresponding to the minimum per 
degree of freedom grow with B in this panel and, for B < 0.1, they fall within the range — 1 = —0.03 ± 0.01. 

In Figures 121 and 131 we plot the adiabatic plus string CMB and matter power spectra superposed using Eq. with 
different values of B and corresponding best fit values of Ug and W for the &-free and 6-fixed models. In both figures, 
the solid line corresponds to i? = and the dash-dot line to = 0,B = 1). 

The simple adiabatic ACDM model with a constant Ug fits the WMAP data very well in general, but much worse so 
on very large scales. For the first few multipoles, 2 < / < 6, the experiment finds a clear deficit of power, as compared 
to the prediction of the model. One could hope that, because cosmic strings generically predict such lack of power 
for low multipoles, adding a string component would improve the fit on large scales. However, as one can see from 
Figs. 121 and 131 this is not the case. The reason is that cosmic strings fit the data extremely poorly in the region of the 
first and the second acoustic peaks, where the WMAP error bars on CiS are very small. The "benefit" from adding 
strings on larger scales, where the error bars are large, is offset by the much larger "damage" that strings cause when 
fitting on smaller scales. 

Clearly, our partial statistical analysis does not yield a rigorous bound on the magnitude of B. However, the results 
in Fig. n can be used to get a rough estimate for the allowed range of the B values. From the discussion at the start 
of this section, a reasonable bound can be based on taking x^(i?)/x^(0) ~ 1-25. In Fig. the maximum value of B 
corresponding to this criterion would be Bmax ~ 0.1. Moreover, we saw that the actual minimum value of x^(^) iii 
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FIG. 2: (Left) Plots of the best fit Cj"^ and ^ computed using tlie &-free model for different values of B together with the 
WMAP's binned data. (Right) Corresponding plots of the galaxy clustering power spectra together with the 2dFGRS data. 
On both plots, the solid line corresponds to i? = 0, dotted line - B — 0.05, short dash line - B — 0.1, long dash line - B — 0.15 
and the dash-dot line corresponds to the pure string contribution, i.e. {W = 0, B = 1). 




FIG. 3: Same as in Fig. |5|but for the 6-fixed model. 



the right panel of Fig.^was at B ~ 0.025, not zero, which indicates that, almost certainly, the available data cannot 
exclude values of S < 0.025. As discussed in S IIII Al this does not simply translate into a constraint on the value of 
the string tension, Gfi. Instead, from Eq. H53|l . B oc (G/i/ A)^. 

Because cosmic strings produce a sizeable vector perturbation, in addition to the _B-polarization, they can also 
induce the B-type polarization of the CMBR. In Fig. 01 we plot the B-polarization angular spectrum, CP^, as 
predicted by our model with B = 0.1. It appears that, even if cosmic strings were contributing only a few percent to 
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FIG. 4: The B-type polarization spectrum C; due to cosmic strings as predicted by our string model with the fit parameter 
B set to 0.1 is plotted as a solid red line. The expected CP^ spectra for E to B lensing (dotted black line) and from primordial 
gravitational waves assuming a tensor-to-scalar ratio of r = 0.1 (dash-dotted magenta line) are plotted for comparison. 

the CMB temperature anisotropy, they could contribute a B-mode signal on scales I ~ 800 that would be in excess of 
the expected B-mode produced by lensing of the E-mode by galaxies |42j| . Such a detection would be an important 
test for the existence of cosmic strings. 



Cosmic strings appear to be a likely by-product of inflation in superstring cosmology. Although these cosmological 
models are only just beginning to be developed in sufficient detail for comparisons with data to be possible, most 
models seem to feature a string scale comparable to or slightly smaller than the GUT scale, and cosmic string tensions 



These expectations motivate studying models in which CMB fluctuations and large scale structure are the conse- 
quence of both growing modes associated with perturbations generated from quantum fluctuations during inflation, 
and active perturbations associated with the gravitational effects of a network of cosmic strings. Here, we have 
presented an attempt at such an analysis, based on a rough comparison of theoretical models with both adiabatic 
perturbations and cosmic strings with WMAP TT and TE power spectra as well as the 2dFGRS results on galaxy 
correlation functions. Our results need to be improved in various ways, both by allowing numerous cosmological 
parameters we held fixed to vary, and by computing likelihood functions, rather than using a rough criterion based on 
a reduced statistic. Nevertheless, they already indicate that although the data currently available generally favor 
models without cosmic strings, they may not exclude nonzero cosmic string tension and density definitively, provided 
the growth of fluctuations in the Universe is not dominated by cosmic strings. Rather conservatively, we have argued 
that values of i? ^ 0.1 are not excluded, although this result needs to be put on firmer footing statistically; we are 
quite confident that B < 0.025 cannot be excluded as yet. Since VB ~ (G.25/A)(G/i/l.l x 10^^), these values of B 
would correspond to 



We have also found the values of rig — 1 that minimze the reduced statistic for different values of B. Conservatively, 
for the range B < 0.1, the various different comparisons with data in the two panels of Fig. ^ are consistent with 
TT-s — 1 ~ —(0.04 — 0.01). Comparing with the predictions of various brane inflation models via Eqs. H22|) and (jSOJ, 
we note that the values we find are in agreement and certainly not in conflict with any of the models. However, we 



IV. DISCUSSION 



Gfi ~ 10-1° - 10-6. 




(57) 
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note that the Coulombic inflation models generally predict |ns — 1| ^ —0.03, at the low end of the range of values we 
infer from comparing with the data, while the powerlaw models allow larger jn,, — 1|; Eq. H22|l implies jn^ — 1| > 0.03 
for any value of ct, and Eq. 1)30(1 implies a constant powerlaw index whose value depends on parameters of the brane 
inflation model. If n<;(fc) can be pinned down with greater precision, it may become possible to discriminate among 
different branc inflation models. We caution, though, since we have not varied the background cosmological model in 
obtaining these bounds, we cannot rule out that a more complete analysis would allow still larger values of G/i, with 
other values of — 1. In a future publication we intend to present results of a more comprehensive study which would 
include varying all cosmological parameters as well as the relevant parameters of the string model, including their 
specific predictions for ns{k), and a better justified statistical analysis based on codes made available by the WMAP 
team for computing likelihood functions for any model for the production of temperature pertubations 0. (Some 
modifications will be needed to account for the nongaussianity of the cosmic string perturbations.) In addition to 
the WMAP and the 2dFGRS data, the analysis will also include the latest results from the Sload Digital Sky Survey 
(SDSS) [13. 

The key to our analysis is the idea that the Universe is a patchwork quilt, with a little bit of cosmic strings 
thrown in to complicate the models. A smoking gun for the existence of cosmic strings could be the detection of B- 
mode polarization on smaller scales at an amplitude considerably larger than is predicted in inflation models without 
cosmic strings. Another potentially distinguishing signature of cosmic strings could be the detection of cosmological 
non-Gaussianity. Tests of the WMAP data have so far been limited to constraining the type of non-Gaussianity 
expected from inflationary models. Some of the most commonly used tests of non-Gaussianity, such as the bispectrum 
test, had actually been shown to be insensitive to possible contributions from cosmic strings [2||. Specially tailored 
tests are likely to be needed to detect string induced non-Gaussianity |43. Cosmic strings may also be detected from 
the observation of identical galaxy pairs in close proximity on the sky j45l |. Gravitational radiation from kinks in 
cosmic strings may also be detectable down to exceptionally small values of G/i ^I^. Pulsar timing is likely to push 
bounds on the density in long wavelength gravitational radiation backgrounds )47j down by an order of magnitude 
or so, corresponding to a factor of three or so in G/Lt/A^ but a substantial contribution to the background from 
supermassive black hole binaries [4^ may frustrate our ability to use these observations to constrain the cosmic string 
tension much further. 
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